Dynamic phase transitions in electromigration-induced step bunching 



o 
o 

(N 
D 
oo 



Vladislav Popko\Q and Joachim Kru^ 
Institut fiir Theoretische Physik, Universitdt zu Koln, Germany. 
(Dated: 6th February 2008) 

Electromigration-induced step bunching in the presence of sublimation or deposition is studied 
theoretically in the attachment-limited regime. We predict a phase transition as a function of the 
relative strength of kinetic asymmetry and step drift. For weak asymmetry the number of steps 
between bunches grows logarithmically with bunch size, whereas for strong asymmetry at most a 
single step crosses between two bunches. In the latter phase the emission and absorption of steps is 
a collective process which sets in only above a critical bunch size and/or step interaction strength. 

PACS numbers; 68.35.-p, 66.30.Qa, 05.70.Np, 81.16.Rf 



o 

CD 



X3 

o 



> 

(N 
(N 
O 

o 



I 

o 
o 



X 



Much of the morphological structure and dynamics of 
crystal surfaces can be understood in terms of the behav- 
ior of steps that separate different exposed atomic layers 
P, 01 • Since they entail a finite free energy cost per unit 
length, steps are long-lived structural defects which nev- 
ertheless, due to their one-dimensional nature, are highly 
sensitive to thermal fluctuations. These fluctuations in- 
duce long-ranged steric interactions between steps, which 
complement similar interactions mediated by bulk elas- 
ticity. When such a system of interacting steps is driven 
out of equilibrium by external forces, e.g. during growth 
or sublimation of the crystal surface, a rich variety of 
morphological patterns and dynamic phenomena emerge. 

As was first shown by Latyshev and coworkers 0, 
step patterns on Si(lll) surfaces can be efficiently ma- 
nipulated by a direct heating current, which induces 
mass transport along the surface through the elec- 
tromigration of adatoms. Subsequently a multitude 
of electromigration-generated step patterns have been 
found and studied experimentally S |^ , including step 
bunches S 1^ , step antibands [3 , in-phase wander- 
ing steps lldl and step pairs many of which still defy 
a comprehensive theoretical description. 

In this Letter we report on the surprising discovery of 
a novel type of phase transition in the most basic model 
of electromigration-induced step bunching originally in- 
troduced by Stoyanov [H E [M Ei . In this model the 
steps are assumed to be straight, and the uniform step 
train is destabilized by an electromigration force in the 
downhill direction. The phase transition occurs as a func- 
tion of a dimensionless parameter &, defined in I^J below, 
which gauges the relative importance of electromigration- 
induced kinetic asymmetry and step drift due to subli- 
mation or growth. This parameter can be tuned exper- 
imentally, e.g., by changing the electromigration force 
through the DC component of the heating current, or 
the sublimation rate through a change in temperature. 

Step drift leads to the exchange of steps between 
bunches, which plays an important role in the evolution 
and coarsening of the bunch pattern 0, . The most 
striking visual signature of the phase transition is a qual- 
itative change in the number of such free steps, and in 
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Figure 1: Typical step configurations (top view of a vicinal 
surface) generated by numerical solution of (0. Graphs (a)- 
(d) correspond to & = 0.1,0.5,5,20 respectively, and in all 
cases U/bl = 0.6. Each frame contains about 60 steps. Re- 
duction of the number of free flowing steps with increasing b 
is evident. In case (d) there are no free-flowing steps since 
the bunch sizes are within the dead zone region of Fig. |3 



the mechanism by which they are exchanged (Fig.^. For 
6 < 1 (strong drift/ weak asymmetry) the step density 
decreases smoothly in the outflow region of a bunch, and 
the number of steps between bunches grows logarithmi- 
cally with the bunch distance. In contrast, for 6 > 1 
(weak drift/strong asymmetry) there is at most a single 
free step between any two bunches, irrespective of their 
size. This feature should make the two regimes clearly 
distinguishable in experiments using reflection electron 
microscopy |3j[3 or scanning tunneling microscopy 6, ^. 

The dynamics in the regime 5 > 1 is remarkably com- 
plex. The exchange of a step is a collective process in- 
volving both the expelling and the receiving bunch, which 
sets in only beyond a critical bunch size, and which is ac- 
companied by breathing oscillations of the entire bunch. 
As a consequence, a stationary bunch shape amenable to 
a continuum description flB | of the type developed pre- 
viously for 6 < 1 does not appear to exist. 
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Model. We consider a system of straight, non- 
transparent steps subject to electromigration and subh- 
mation (including also a growth flux is straightforward) . 
We work in the attachment-detachment limited regime, 
where the kinetic length d = D/k, the ratio of surface 
diffusion coefficient D and attachment rate k, is large 
compared to the step spacing 0. The equations of mo- 
tion for the step positions Xi{t) then take the form 



dxi 
~dt 



Xi) + 



+C/(2/,-/,_i 



(1) 



where the time scale has been normalized to the sublima- 
tion flux. Summing over i we see that the average step 
velocity v is equal to the mean terrace width I. In numer- 
ical solutions of lengths are measured in units oil, i.e. 
we set 1 = 1. The last term on the right hand side rep- 
resents stabilizing step-step interactions of strength C/, 
where, for combined entropic or dipolar elastic repulsion. 
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with V = 2 1]. The parameter b governs the asymmetry 
between ascending and descending steps, relative to the 
mean step velocity, which induces step bunching when 
6 > 0. Linear stability analysis of Q shows that the in- 
stability sets in at wavelengths corresponding to bunches 
containing more than M* steps, with 



M* = 27r[arccos(l - bl/l2U)Y 



(3) 



In previous work more complicated variants of the step 
equations have been studied numerically, and some of 
the features analyzed in this paper have been described 
on a qualitative level 0, 0| . The advantage of using 
the attachment/detachment limited dynamics (^3) lies in 
the linearity of the destabilizing terms, which allows to 
clearly expose the key role of the parameter b and the 
existence of a sharp phase transition. 

In terms of physical quantities, the parameters b and 
U are given by |13| 
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where T is the step mobility for the Brownian motion of 
an isolated step is the atomic area, F is the elec- 

tromigration force acting on an adatom, Tg is the inverse 
desorption rate, a = a/l is the miscut angle, and g is the 
step interaction parameter. 

The model Q is expected to apply in two of the four 
temperature regimes |j] in which step bunching is ob- 
served on Si(lll), around 900° C and around 1250° C 



0, ■ The parameters given in |^ lead to the esti- 
mates & « 14 in the low temperature regime and h k, 0.3 
in the high temperature regime, which shows that both 
cases 6 < 1 and 6 > 1 are experimentally realizable. 

Step equations of the form can also be derived for 
step bunching induced by Ehrlich-Schwoebel (ES) bar- 
riers during sublimation |l7l | or by inverse ES barriers 
during growth 0, 0|. In this sense ^ constitutes a 
rather generic model of step bunching kinetics. However, 
in step bunching induced by ES barriers the parameter 
b is restricted to the interval < 6 < 1, and hence the 
phenomena described in this paper do not occur. 

Structure of the outflow region. In the presence of step 
drift, coarsening of step bunches is a very dynamic pro- 
cess during which steps continuously leave (flow out of) 
one bunch and join (flow into) its neighbour 

[ilEi In 

[lil Il9| it was shown that the analysis of the outflow re- 
gion provides key insights into the shape and dynamics 
of bunches for 6 < 1. We shall see now that there are 
drastic differences between the outflow regions for the 
cases & < 1 and 6 > 1. We consider a bunch contain- 
ing a large number M 1 of steps, so that its shape 
can be considered quasi-stationary. We impose periodic 
boundary conditions Ai(t) = Ai_|_M(i) for the terrace 
sizes A; = Xi+i — Xi. Stationarity implies then periodic- 
ity of each step trajectory (up to an overall shift with ve- 
locity V = I), with some period r(&, U,M), during which 
each step i will once cross the plateau between two con- 
secutive bunches. After time t/M, each step i will sub- 
stitute the position of step i + 1 (up to a constant shift 
independent of z), so that 
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Deriving an equation for A^ (t) from ^ and substituting 
we get a differential-difference equation for a single, 
T-periodic function A{t) = Ai(t) 18] 



dA{t) 
dt 



-&A(i)- 



-Ui. 



(6) 

where for brevity the J7-containing terms are only 
sketched. The (unknown) period r determines the veloc- 
ity of a bunch: after time r the bunch shifts by (—Ml) 
in a frame co-moving with velocity v; in the laboratory 
frame its lateral velocity is then 



V = l{l- M/t). 



(7) 



Big bunches are separated by wide plateaux, and for the 
steps crossing a plateau (in case there are many) the U- 
term in lO should become negligible. In this outflow 
region one can solve the remaining linear part by the 
ansatz A{t) ~ exp^qMt/r) obtaining the transcendental 
equation 6 (cosh (5) — 1) = sinh (g) — qM/r. To fix the 
unknown parameter r, we recall the Fourier analysis of 
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Figure 2: Full line: Distance A(t) between a pair of neighbor- 
ing steps moving through the bunch interior (initial and end 
regions excluded), for b = 10, U/l = 6, M = 64. Symbols show 
the sizes of 42 consecutive terraces (out of 64) , at times t — 
(+), t = 0.34 • (r/M) (x), t « 0.55 • (r/M) (*), illustrating 
the oscillatory breathing of the bunch. They lie on the curve 
A{t) because of ©. Inset: A(t) for b = 0.176, U/l = 0.108. 
Oscillations are triggered by a step colliding with the front 
end of the bunch, but do not extend into the outflow region. 



Figure 3: Phase diagram characterising quasi-stationary step 
bunches of size M as function of step-step repulsion U for 
6 = 10, / = 1. Three phases can be defined: (a) Bunches 
smaller than the size Ai* deduced from linear stability anal- 
ysis [Eq.(|^] dissolve (filled circles). (b) Dead zone: for 
Mci{U) < M < Mc2{U) bunches are stable but do not emit 
steps (squares), (c) For M > Mc2{U) and M* < M < 
Mc\{U) bunches are stable and emit steps (crosses). The lines 
Mc\{U),Mc2{U) terminate at a critical interaction strength 
Uc{b) beyond which stable bunches always emit steps. 



(Q) in which shows, irrespective of the value of 6, 
that for large M generically t{M) « M -f 0(1). Thus 

6 (cosh (7 — 1) = sinh(7 — q (8) 

which has a real positive solution for each b < \ but 
no solutions for 6 > 1. In the following we explore the 
consequences of this fact. 

Number of steps between bunches. For 6 < 1, the 
existence of a solution (7 of implies a smooth de- 
crease of the step density in the outflow region, with the 
terrace widths increasing exponentially, as Ak/Ak-i ~ 
exp {q). To estimate the number Nf of free steps between 
two bunches of size M, we equate the total length ~ 
I exp[giV/] occupied by Nf terraces to the typical distance 
Ml between the bunches, and obtain Nf « q~^\nM. 
For small b, the solution of © can be approximated by 
q « 35. 

For b ^ 1, q diverges and Nf vanishes. The absence of 
solutions of (IHJ) means that the J7-term in can never 
be neglected and that correspondingly there can be at 
most one step crossing the plateau between two bunches, 
at any stage of evolution. One can check, using ^ and 
(O, that any configuration with more than one step be- 
tween two bunches is unstable for 6 > 1, so that all steps 
except at most one will be pushed back to the bunch they 
originated from. In Fig. ^ we show numerically generated 
bunch configurations in the course of coarsening for 6 < 1 
and 6 > 1, which confirm this conclusion. 

Dynamics of emission and absorption of steps. We 
now examine in more detail how steps are emitted from 
a bunch. It is seen directly from that the last step 



(with label i, say) of the bunch at position Xi, which is 
trailing a wide terrace of width = x^+i — l^ will 
be driven to escape from the bunch by the linear term 
(1 — 6)Ai/2 , provided 6 < 1. This term indeed gives the 
main contribution to the dynamics of the last step of the 
bunch, as we see from numerical analysis. The emitted 
step does not perturb the remaining steps; the (i — l)-th 
step which has become the last, is free to escape once 
the i-th step has travelled sufficiently far. Bunches emit 
steps continuously, creating an outflow region governed 
entirely by the linear part of (^) 0, 0| . 

In contrast, for > 1, the linear term (1 — b)Ai/2 
in ^ gives a negative contribution to the step velocity, 
and the only way to move the last step i away from the 
bunch is by step-step interactions [the [/-term in (^]. 
Since the next step i — 1 cannot be emitted before step i 
has landed at the next bunch, the conflguration of steps 
at the end of the bunch has to be changed by the emission 
process - if it were unchanged, the next step i — 1 would 
be emitted immediately after the i-th. This gives rise to 
oscillations of the bunch proflle at the end, which spread 
to the whole bunch, and whose amplitude grows with 
increasing b. Such oscillations at the outflow end of the 
bunch are completely absent in the 6 < 1 phase (Fig. 

When the emitted step flnally collides with the receiv- 
ing bunch it provokes perturbations in the inflow region 
of the bunch, which are visible both for 6 < 1 and b > 1. 
In the case & < 1, however, the oscillations in the in- 
flow region are damped and disappear towards the inte- 
rior of the bunch. On the contrary, in the b > 1 phase 



4 



100 



90 




10 20 30 40 50 60 70 
t 



Figure 4: Space-time plot of individual step trajectories dur- 
ing coarsening, in a frame co- moving with the step velocity v, 
for b = 20, U/l — 12. Step emission starts when the bunch size 
exceeds a critical value. The isolated emission event around 
t ~ 7 occurs because the bunch shapes are not yet station- 
ary. The initial configuration is a slightly perturbed train of 
bunches consisting of 16 steps each. 



the oscillations penetrate through the bunch, regain their 
large amplitude towards the bunch tail, and culminate in 
the emission of the last step of the bunch, provided that 
the initial impact was sufficiently strong (Fig. [5J • The 
persistence of oscillations through the bunch interior im- 
plies correlations between the emission and absorption 
of steps, which should have important consequences for 
the coarsening dynamics; this question will be addressed 
elsewhere. 

Onset of step emission. We have seen above that the 
emisssion of steps in the 6 > 1 regime is a nontrivial 
dynamical phenomenon facilitated by a large step-step 
repulsion U, and suppressed by the kinetic asymmetry b. 
For small U (or large b) the oscillatory breathing of the 
bunch may not be able to trigger the emission of steps 
when bunches are small. The typical behavior of bunches 
as a function of size M and step interaction [/ at a fixed 
value of 6 > 1 is summarized in the phase diagram in 
Fig. 13 For any given b > 1, there exists a critical value 
Uc{b) such that for U < Uc bunches emit steps only for 
sizes M* < M < Md and M > Mc2, whereas for U > Uc 
stable bunches always emit steps. Inside the dead zone 
Mci < M < Mc2 the time interval t/M between emission 
of steps is infinite, and correspondingly bunches move 
with the mean step speed, V = v [see (OJ. The ratio 
t/M decays monotonically to 1 with distance from the 
dead zone, reflecting the fact that limM^cxj t/M ~ 1 for 
any fixed b 0. 

Diagrams for different b can be superimposed after 
rescaling Uc and Ale = Mci,2{Uc) according to the rela- 
tions Uc « 0.0105 • 6" and Mc ~ 2.112 • b'' with a « 2.87, 
7 w 0.935 for all parameters investigated (3 < 6 < 25, 



10~^ < U < 120), with a relative error not exceeding 5%. 
Note that the relation 7 = (a — 1)/2 implies invariance of 
the linear instability curve ((SJ at large U under rescaling. 

Different step kinetics for bunches of different sizes 
implies a change in coarsening dynamics, highlighted in 
Fig. 0] For b:^ 1, depending on the value of U different 
coarsening scenarios are possible. For U > Uc steps are 
exchanged throughout the coarsening process, while for 
U < Uc/2 late stage coarsening proceeds in two stages: 
without step exchange (for bunches sizes smaller than 
Mc2) and with step emission once the typical bunch size 
exceeds Mc2- Coarsening with or without step exchange 
has previously been associated with nonconserved (b fi- 
nite) and conserved (6 = 00, no sublimation) dynamics, 
respectively ; here we see that both types of behavior 
may coexist when b > 1. 

Conclusions. We predict a new type of phase tran- 
sition in electromigration-induced step bunching within 
the regime of nontransparent steps and attachment- 
detachment limited kinetics. The transition is charac- 
terized by a dramatic change in the number and behav- 
ior of the free steps that are exchanged between bunches, 
which should be clearly visible in experiments on surfaces 
vicinal to Si(lll). Theoretical challenges for the future 
include the development of a continuum description for 
6 > 1, and the investigation of the correlated coarsening 
dynamics in this regime. 
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KR 1123/1-2. 



* E-mail: popkov@thp.uni-koeln.de 
^ E-mail: krug@thp.uni-koeln.de 
[1] H.C. Jeong and E.D. Wilhams, Surf. Sci. Rep. 34, 171 
(1999). 

[2] J. Krug, Int. Ser. Num. Math. 149, 69 (2005). 

[3] A.V. Latyshev, A.L. Aseev, A.B. Krasilnikov, and 

S.I. Stenin, Surf. Sci. 213, 157 (1989). 
[4] K. Yagi, H. Minoda, and M. Degawa, Surf. Sci. Rep. 43, 

45 (2001). 

[5] H. Minoda, J. Phys.-Cond. Matter 15, S3255 (2003). 
[6] Y.-N. Yang, E.S. Fu, and E.D. Wilhams, Surf. Sci. 356, 
101 (1996). 

[7] J.J. Metois and S. Stoyanov, Surf. Sci. 440, 407 (1999). 
[8] K. Fujita, M. Ichikawa, and S.S. Stoyanov, Phys. Rev. B 

60, 16006 (1999). 
[9] K. Thurmer, D.-J. Liu, E.D. Wilhams, and J.D. Weeks, 
Phys. Rev. Lett. 83, 5531 (1999). 
[10] M. Degawa, H. Nishimura, Y. Tanishiro, H. Minoda, and 

K. Yagi, Jpn. J. Appl. Phys. 38, L308 (1999). 
[11] O. Pierre-Louis and J. -J. Metois, Phys. Rev. Lett. 93, 

165901 (2004). 
[12] S. Stoyanov, Jap. J. Appl. Phys. 30, 1 (1991). 
[13] D.-J. Liu and J.D. Weeks, Phys. Rev. B 57, 14891 (1998). 
[14] S. Stoyanov and V. Tonchev, Phys. Rev. B 58, 1590 
(1998). 

[15] M. Sato and M. Uwaha, Surf. Sci. 442, 318 (1999). 
[16] A. Pimpinelli, V. Tonchev, A. Videcoq, and 
M. Vladimirova, Phys. Rev. Lett. 88, 206103 (2002). 



5 



[17] J. Kmg, V. Tonchcv, S. Stoyanov, and A. Pimpinelli, [19] F. Slanina, J. Krug, and M. Kotrla, Phys. Rev. E 71, 

Phys. Rev. B 71, 045412 (2005). 041605 (2005). 

[18] V. Popkov and J. Krug, Europhys Lett. 72, 1025 (2005). 



